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Pe oka T 


The problem of extending vertical tidal and river datums 
through the ocean/river transition zone of navigable rivers 
is examined through the analysis of water-level measurements 
from an ocean anda river station at San Francisco and 
Sacramento, California. Modification of the tides after 
passage into the transition zone was analyzed by decomposing 
the raw water-level data for both stations into a tidal and 
a non-tidal component, and comparing the tidal components 
for tide range ratios and time differences. It was deter- 
mined that as the mean river stage increases, the range ratio 
and the effective tide wave speed both decrease and the sym- 
metry of the tide wave changes. Of the six standard Pacific 
Coast tidal datums and five river datums defined in this 
study, MSL-MRL is the only common tidal/river datum that is 
continuous through the transition zone. The MLLW ocean chart- 
ing datum and the MHW tidal waterfront property boundary datum 
can be carried upriver, where both effectively merge with the 


MRL, by a separation-addition procedure. 
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CS ROOUCT LON 


Vertical datum planes are used routinely for charting 
purposes and for the determination of waterfront property 
boundaries. In the fully tidal marine environment datum 
peeicoeeare well-defined and procedures for their determina- 
tion are well-established (Marmer, 1951; Shalowitz, 1964). 
As these planes are extended up a navigable river they be- 
come less well-defined. Within a short distance the magni- 
tude of the annual river stage variation becomes the same 
order as the range of the tide; further up the river the 
river stage variation may be much larger than the tidal 
range. Unlike the tides, which are highly repetitive and 
lend themselves to accurate datum determination, river levels 
are highly irregular in both height and duration and require 
meet rerent approach for datum determination. Accordingly, 
datum planes that are useful in the purely tidal environment 
are not necessarily useful or desirable in rivers, and vice 
versa. The primary purpose of this study is to explore the 
usefulness of vertical datum planes in the transition zone 
between the fully tidal marine regime and the non-tidal 
river regime. 

In order to understand what happens to conventional 
tidal datum planes as they are extended up a river it was 
necessary to determine how the tides change in response to 


the river stage variations. In order to accomplish this, 





two water-level measuring stations on the same river system 
were chosen. The first was a fully tidal station, the 
National Ocean ames reference tide station at the Presidio, 
Poeeraneisco, California. The second was the U.S. Geologi- 
cal Survey gaging station at Sacramento, California in the 
tidal and riverine transition zone. Both stations may be 
considered to lie on the Sacramento River, as shown in Fig- 
ure 1. A 19-year hourly water-level time-series record for 
the tidal epoch 1959-1977 was obtained for each station and 
analyzed. 

The river stages at Sacramento vary from heights of 
above 20 feet on the water-level gage during ordinary winter 
and spring runoff to less than 2 feet during the dry summer 
season. Visual observations of the raw water-level record 
Mietedre that the tidal influence there is greatly diminished 
or absent during high river stages. In order to examine the 
character of the tides at Sacramento, the hourly raw data 
were decomposed into a tidal and a non-tidal component; this 
was done for San Francisco as well. Figures 2 and 3 show a 
Weaeeor raw water-level data for both stations, the tidal 
component extracted from the same data, and the residual 
water level. The range of the tide at San Francisco and 
Sacramento obtained from the tidal component was then com- 
pared and the time difference between the tide wave arrival 
at the two stations was determined,both as a function of the 
river stage at Sacramento. 

The Sacramento water level remaining after the tidal 


component is removed is considered here to be representative 
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of the pure river regime found upriver from Sacramento at 
Beene remasbeyord the effective reach of the tides, whereas 
the raw water-level record occurring at times of succes- 
sively lower river stage 1S considered to reasonably repre- 
sent the successively increasing oceanic environment down 
river from Sacramento. Thus, in effect, all degrees of 
feast elon from full riverine to full tidal environment 
are represented in the water-level data from the two stations. 
The results of this investigation are intended to provide 
technical guidance for the definition of datum planes in the 
ocean-river transition zone. Emphasis is placed on datums 
used for navigation and property boundary determination. [It 
does not consider the legal aspects of property boundaries 


mmene transition zone. 
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ii. DATA ANALYSIS AND PROCEDURES 


The raw data for this investigation were supplied by the 
National Ocean Survey and the U.S. Geological Survey. They 
consisted of a 47-year time-series of hourly water-level 
measurements at the Presidio, San Francisco; and a 16-year 
time-series of 15 minute water-level measurements and a 4-year 
time-series of hourly water-level measurements at the U.S. 
Geological Survey gaging station at Sacramento (Table 1). 

The San Francisco data are referenced to a staff zero which 
is 8.61 feet below the National Geodetic Vertical Datum 
(NGVD), formerly termed Sea Level Datum 1929 CSLD 1929). The 
sacramento data are referenced to a staff zero at 0.00 feet 
NGVD COltmann, 1378, personal communication). For consis- 
Geney in handling, all data were converted to a digital in- 
Memval Of One hour. The hourly data are available at the 
Department of Oceanography, Naval Postgraduate School on 
computer cards and at the Oceanographic Division of the 
National Ocean Survey on magnetic tape. 

In order to understand the analyses performed, some kKnow- 
ledge of the tide and river-level characteristics at both 
Pome cencisco and Sacramento is desirable. The tide at San 
Prancisco is a mixed tide. The range of the tide varies from 
about five feet during spring tides to about three feet during 


neap tides. The elevations of the standard tidal datum planes 


ie 





and the geodetic datum plane at San Francisco for Epoch 1941- 


fo9 are; 
Ref to Ref to Ref to 
Start 70 NGVD MLLW 
MHHW id tema ec ahs or ae See le gam 
MHW UO ee sais: Lai. ae ae 
MTL Seco O25 Sle 
MSL 8.80 — 0.19 Be05 
NGVD como E Cow 2ue0 
MLW oy -1.74 ea? 
MLLW Oras —Vigenells Ce 010 
The tide at Sacramento is also a mixed tide. Its range 


at spring tides varies from about three feet at low river 
Stages to under one-tenth of a foot at high river stages. The 
river level, after the tides have been removed, is highly 
variable, with extreme values ranging from under two feet to 
greater than 29 feet during the period studied (1953-1977). 
Appendix B contains graphical printouts of these data. 

In order to understand the behavior of datum planes that 
extend into the transition zone from both the fully tidal and 
fully non-tidal regimes, it was necessary to analyze the 
water-level data to determine the modification of the tides 
after passing into the riverine environment. The analyses 
performed of the hourly tide data at both San Francisco and 
Sacramento involved the following sequential steps: (1) De- 
composition of the raw hourly time-series data into a tidal 


constituent and residual water levels in time-series form; 


Ihe 








(2) determination of the height and time of each tide in the 
extracted tide curve by means of a fitted polynomial eo Ene 
high and low Teas - (3) comparison of the tides at each sta- 
tion to determine the effect of river stage on the heights 
and times of high and low waters; and (4) analysis of the raw 
time-series record at Sacramento for use in determining river 
datum planes. The non-tidal component at Sacramento was 

also analyzed and compared with ee ares in (4) above. 


A description of each analysis and the reason for it is given 


below. 


peel eCOMPOSITION OF WATER-LEVEL RECORDS 

In order to compare the ranges and times of the tide at 
the two stations, the raw data were decomposed into a tidal 
and a non-tidal component. The primary tidal constituent has 
a period of 24.8 hours; accordingly, 1t was decided to compute 
a 25-hour running mean to accomplish the primary separation. 
Figures 4 and ° show an example of the raw and decomposed 
data at San Francisco and Sacramento, respectively, for a 
given five-day period. The computer program for this separa- 
tion is included in Appendix A-1l. This data reduction and 
all other computational work was accomplished on the IBM-360 
Pemetner act the W. R. Church Computer Center of the Naval 
Besepraduate School. 

This separation procedure proved to be effective in fil- 
tering out regular variations with a period of 25 hours or 
less. Two distinct cyclical components remain in the Sacra- 


feomeemtime-series. The dominant One 15 the irregular annual 
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variation in the river level due to seasonal precipitation 


moceoumort. Inhis variaticn is of a magnitude Ornup orc! 
feet per year. The other is a harmonic component with a 
fortnightly period which appears to be linked to the spring- 
Meap tide cycle. This component has a magnitude of less 
than 0.4 foot and was not filtered from the Sacramento re- 
Sidual river level. McDowell and O'Connor (1977) describe 

a rise in mean river level for wean movers ian Lndia 
caused by an accumulation of water upriver during the spring 
tides, and evidently it is this phenomenon that remains in 
the Sacramento time-series record. Pigure 6 displays this 
meeumientliy cycle. The cycle evidently does not occur in 
the essentially fully marine environment at San Francisco. 

It may be noted from Figure 5 that the tide waves at 
sacramento are markedly asymmetrical, and in striking con- 
trast to the tides at San Francisco for the same period 
shown in Figure 4. The asymmetry is attributed by McDowell 
emewe Connor (1977) to the fact that the tide crest moves 
faster than the tide trough due to increased water depth at 
high tide compared to that at low tide. The high tide marked 
feeis a Common tide in both figures. 

In the San Francisco time series, the water-level varia- 
tions that remain after the removal of the diurnal and 
shorter tidal components are irregular and of very low ampli- 
tude. These variations are attributed principally to meteo- 


rological effects, but must also contain very low amplitude 
long period tidal constituents, although none are visually 


evident in graphical printouts of the residuals (Appendix B). 
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B. HIGH WATER AND LOW WATER HEIGHTS AND TIMES 


Comparison of tide heights at San Francisco and Sacra- 
mento cannot be done meaningfully in the usual way as with 
ocean stations. The raw water-level data at a river station 
Semedins not only the tidal component, but also the river 
level, which may have any arbitrary height. In order to 
make a comparison, the high/low water heights and times 
needed to be extracted from the tidal component time-series 
record. Although heights can be obtained accurately direct- 
ly from hourly digital data by use of a computer, determina- 
tion of the times of high/low waters needed to be refined in 
order to obtain the time to an accuracy better than the 
elosest hour. This was accomplished by curve-fitting to the 
data points to obtain the times and heights to the desired 
precision of 0.1 hour and Q.1 foot. There were also many 
irregularities in the time-series record not related to the 
tide, especially in the Sacramento data, which made curve- 
meerame additionally desirable. 

Curve-fitting was accomplished using a computer program 


Mmeaberen for this investigation which is best understood by 


reference to Figure 7. The first step was to determine an 
Mmemcial suess as to the time of the high or low water. The 


tide heights at hour times t, and t, were summed, then sub- 
metetodetrom the sum of the heights at times tf, and t.. The 
algebraic sign of this difference was noted and associated 
With the mid-point, time t- This procedure was indexed for- 


ward one hour and repeated, comparing heights at times t., 





emcee With the heights at times Pec tne Aleem ralc 
sign of the difference was associated with time the This 
procedure was Laclos eel an hour at a time throughout the entire 
record. Whenever the algebraic sign changed it was considered 
that a high or low water was associated with the mid-point 

Mee nourly data point) of the last summation. This initial 
guess procedure gives the time of the high or low water to 

mae eCllosest hourly data point. er sae eee A-2 contains the 
Semputer program. 

Once an initial guess was made for the high/low water 
time, the associated water level and the two preceding and 
following hourly water levels were fitted to a second-order 
polynomial. The fitting was done in a least squares sense. 
Since the data points are at a regular one-hour interval, 
by adjusting the indexing the normal equations for the least 
Squares fit needed to be solved only once. The heights of 
the five data points are simply substituted into a single 
equation for each of the three coefficients of the quadratic 


equation: 


2 


‘Ge= 2 x 


+ a.X +f 4 
¢] ils 


2 


where Y is the height of the actual high/low water and X is 
miemedme Of the actual high/low water. This quadratic was 
Prema tirepoentiated, set equal to zero, and solved for 4, 
epemactial time Of the high/low water. This time Cto the 
plesesteasl hour) was substituted into the quadratic to obtain 


the height of the high/low water to 0.1 foot. By way of 
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example, Figure 8 shows the tide heights and times derived 
from this procedure superimposed upon the raw data water- 
levels for San Brariedsc6 and Sacramento (shown in Figures 4 
and 5). Table 2 compares the heights and times of the high/ 
Mow waters tor the San Francisco station calculated by this 
procedure with the heights and times derived manually by NOS 
from an analog tide record for the same dates. Appendix A-3 
Semrains the normal equations fom ORE least squares fit and 
aiso the computer program. 

(ies op ]ective curve-fitting technique has two potential 
inadequacies which are dependent upon the raw data. First, 
Mempeme Paw data are not sufficiently smooth, extra “high/low 
waters” may be extracted which represent small irregularities 
in the data. Each high or low series of points in the data 
which persists for two or more hours produces a high or low 
water. The second inadequacy occurs when the range of the 
Beeeettde 1s sufficiently small, in which case the high or 


low waters go undetected. Both of these cases occurred in 
the Sacramento data, and only the first case was observed in 
the San Francisco data, but only very occasionally. The 
handling of these problems is described below. 

The heights and times of the high and low waters at San 
Francisco and Sacramento were used to determine how the tide 
varies in response to the river stage at Sacramento. Other 
river factors, such as time-rate-of-change of mean river 
height, were considered to be of second-order importance and 


were not examined. 


ims 





The range of the tide rather than the heights of the 
high/low waters at the two stations was chosen for comparison. 
The latter procedure is conventional when the stations are 
both fully tidal, but ceases to be satisfactory where both 
the height and time of the tide vary significantly in re- 
sponse to a varying river height. Accordingly, the range of 
the tide at Sacramento was divided by the range of the tide 


at San Francisco for the same tide,.to give a range ratio, 
RO 
iimeerder (oO properly match up the tides at both stations, 
and in view of the fact that a significant number of the 
tides were missed and some false tides were recorded at 
sacramento, the high and low waters were filtered to ensure 
that only genuine tides were compared. The first filtering 
utilized the time difference between the tide-wave passage at 
San Francisco and its passage at Sacramento. Tide tables and 
hand analysis of the water-level data indicated that this 
@mererence is about nine hours. After initially matching the 
tide by hand, the computer program would look at the time lag. 
[ima ende at Sacramento occurred less than 6.5 hours after 
the San Francisco tide, it was assumed that a fictitious 
tide had been detected at Sacramento. The Sacramento tide 
Peeecteregarded and the next Sacramento tide compared. If 
the initial time difference was greater than 14.5 hours it 
Meeeassumed that a tide had been missed at Sacramento, so the 
San Franciseo tide was disregarded and the next San Francisco 


[id -meompared. The 6.5 and 14.5 hour cutoff points were arbi- 


trarily chosen after hand analysis of the tide records. 


ike 








The second filter developed was required because of the 
meeoeantial number of real and false ides with a range mag- 
nitude on the order of one-tenth of a foot or less at Sacra- 
mento. These were especially prevalent when the river level 
was greater than 1l feet. This second filter disregarded 
all tides with a range less than 0.1 foot. This difference 
Beeeeen= Same order of Magnitude as "noise™ in the data. 

A third filter was used to eliminate the few incorrect 
Pange comparisons that survived the first two filters. Since 
— contains an algebraic sign as well as a magnitude, if the 
range ratio was not a positive value, a falling tide was ob- 
viously being compared with a rising tide, or vice versa. 


Any Ae values that were negative were discarded. 


GeeeexANGE RATIOS AND PHASE INTERVALS 

ee thes point the individual R_ values retained were con- 
sidered to be valid. Visual examination and hand analysis of 
mgemdaca determined that both R. elide thew time difference, “At, 
between a tide passage at San Francisco and its passage at 
Sacramento are dependent on the tide phase being considered. 
That 1s, the average value of R, for the tide phase LLW-LAW 
is different from the average value of RS for the phase LLW- 
HHW for the same river level. Also, At for a LLW is different 
from At for a LHW for the same river level. Accordingly, in 
order to determine which of eight possible tide phase inter- 
vals a measured 8, is associated with, or which of the four 
possible tide phases a At is associated with, at least five 


consecutive tides covering four successive phase intervals 
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must be detected by the computer program. This allows deter- 
mination of the tide type (HHW, LHW, etc.) and the phase in- 
terval (CHHW-LLW, HHW-HLW, ete.) for the river level occurring. 
If five consecutive tides could not be compared, tides were 
discarded until five consecutive tides were available. Since 
Smpiase Change occurs about the time of neap tide, one or two 
R, values about the time of neap tide are discarded by this 


procedure. This procedure is applied after all of the other 


filters. 

Nine years of overlapping water-level data were analyzed 
in this manner and approximately 34% of the possible common 
tides at San Francisco and Sacramento survived the above fil- 
ters. This amounts of 4392 data points. The average range 
ratios, RO and time differences, At, for all phase intervals 
mee eaae types are displayed in Figures 9-11 and Tables 3 


eu 4 . 


D. MEAN WATER-LEVEL DETERMINATION 

The Sacramento time series was analyzed both by calendar 
years and water years (1 October-30 September). Water years 
are commonly used in river studies in the United States in 
order to avoid splitting the peak precipitation season (winter) 
between two years. Calendar years are commonly used in tidal 
studies. Five river datum planes were defined and calculated 
mormeene 19-year period 1959-1977. They include the Mean 
River Level (MRL), the Mean Flood Level (MFL), the Mean Low 
River Level (MLRL), the Mean High Monthly River Level (MHMRL), 


and the Mean Low Monthly River Level (MLMRL). 
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The MRL is the 19-year mean of the average water level 
memeecach individual calendar or water year determined from 
raw hourly heights, each year being given equal weight. Each 
year represents 8760 hourly data points (8784 during leap 
years). The MFL is the 19-year mean of the single extreme 
high raw hourly water level for each year. The MLRL is the 
19-year mean of the single extreme low raw hourly water level 
for each year. The MHMWL is the 19-year mean of the Single 
highest month in each year determined from raw hourly heights; 
the highest month is not always the same calendar month in 
each year. Similarly, the MLMRL is the 19-year mean of the 
Single lowest month in each year. The latter two datums, 
chosen arbitrarily, were defined because they represent 
moderately high and low river levels lying between the MRL 
and the extreme high and low river datums, MFL and MLRL. It 
eijeula be noted that other river datums could be defined. 

MemeeGaln an indication of the vertical variability, or 
meee lity, of these planes, a 95% confidence interval was 
Meeermined for each datum using the Student-T distribution. 
With a confidence of 95% the true mean lies within plus or 
minus one interval of the calculated mean. The standard de- 
mee toms Lor each of the datum planes were also calculated. 
The mean yearly and 19-year elevations of these five river 
ieaims. alone with the deviation of the yearly values from 
the 19-year mean, are shown in Figures 12-15. The 19-year 
Peeawere also given in Table 5; the tidal datums at San 


Francisco are included for comparison. 
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A separation-addition procedure for constructing river 
datums, which is explained in a later section, requires de- 
Semposition of sere water-level data from Sacramento into 
a tidal and a non-tidal or river component, followed by com- 
putation of the mean river level derived from either the raw 
or the non-tidal hourly water levels. Accordingly, compari- 
son of the mean monthly and yearly water levels at Sacra- 


mento for the period 1958-1963 was made using both the raw 


hourly river data and the non-tidal component. The results 
are given in Appendix C. The comparative values are seen to 
@eeee to Within 0.02 feet in most cases. Thus, to a high 


degree of accuracy, those river datums that are derived from 


hourly readings may be obtained from either the raw data or 


from the non-tidal river component. 
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ite eer rECTSOF RIVER STAGE UPON 
PEPER CiaRAC TEREST ics 


It has been determined from nine years of comparative 
Mepeer fevels (1952-1961 and 1967-1971) that, in general, as 
the river stage at Sacramento increases, the range ratio be- 
tween Sacramento and San Francisco decreases and the time 
@meterence increases. These results are shown in Figures 9- 
ll and in Tables 3 and 4. It may be seen that the range 
Maelo decreases approximately exponentially. The time dif- 
ference increases approximately linearly at low river stages 
and tapers off at higher river stages. The reader should 
note that the numbers of values in Table 4 at river levels 
greater than or equal to ten feet are too small to give re- 
Meesemtative plots of At in Figure 11. River level, as used 
in this study, refers to the residual river level remaining 
maeeeeeemOvVal Of the diurnal tides. 

Two thousand two hundered cases of falling tide were ex- 
amined. As may be seen in Figure 9 the relationship between 
RL e@ammsiver level is very similar for all four phase types. 
The mean value of oe varied from 0.249-0.292 at river levels 
momeenie range of 2-4 feet to 0.026-0.034 at river levels of 
met feet. A very small number of river levels Delow two 
feet were observed and they are included in the 2-4 foot cate- 
gory. Two thousand one hundred ninety two cases of rising 
tide were examined (rigure 10). — varied irom 60212-07331 


at river levels of 2-4 feet to 0.026-0.034 at river levels or 
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iz—is teet. As stated previously, staff zero at Sacramento 


lies at 0.00 feet NGVD. 

ivyemcaousand he hundred cases of low waters were ex- 
amined (Figure 11). The mean time difference of all low 
waters was found to increase from 9.14 hours at river levels 
in the range of 2-4 feet to 9.84 hours at river levels of 
12-14 feet. Two thousand one hundred ninety two cases of 
high waters were examined. The teem time difference in- 
creased from 8.63 hours at river levels of 2-4 feet to 9.31 
hours at river levels of 12-14 feet. The rate of increase 
in At for both low and high waters 1s approximately ten 
minutes per foot increase in river level up to a river level 
of 7 feet, and substantially smaller at higher levels. It 
may be noted that at low river stages the NOS Tide Tables 
Maeaeate a difference of 7.30 hours for high tides and 9.28 
hours for low tides. This is in close agreement with this 
eemcdy for low tides, but differs by approximately one hour 


fem hreh tides. 
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> iii VERTICAL DATUMS 


For a vertical datum plane to be useful it should possess, 
Memeat=s SPeat a dégree as possible, the following properties: 
(1) maximum vertical stability, which requires water-level 
measurements over as long a period as possible; (2) easy visu- 
alization, which requires a oe definition; (€3) derived 
from continuous or suitably digitized time-series measure- 
ments at a water-level gage; (4) obtainable from simple sta- 
feotteal handling procedures and not from theoretical or 
numerical models; (5) the datum should run through the ocean/ 
river transition zone without discontinuity; and (6) the 
datum plane, within the riverine environment, should lie 
within the normal channel of the river. The two most impor- 
meatier Or these properties, stability and continuity, are 
discussed below. 

With regard to identifying datums that can be carried 
mico and through the transition zone, the standard oceanic 
feeum planes derived by NOS from the San Francisco tide 
measurements were considered: MHHW, MHW, MTL, MSL, MLW, and 
MLLW. For the river, five datums derived from river stage 


measurements defined previously were calculated: MFL, MHMRL, 


rt 


Fy 


MeeeeMiMei, and MLRL. The river planes were all computed fo 
maeeus=vear interval 1959-1977. Nineteen years was chosen 


Pempeme@eonsictent with a tidal epoch, although there are no 


variations apparent in the river levels of that periodicity. 
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itesriver planes are shown in Figures 12-15 and will be 
@ensrdered first in regard to vertical stability. In Table 
5 their stability is given in the confidence intervals and 
standard Pevaattons. solely from the standpoint of stability, 
the two low water datums would be the most desirable. How- 
ever, the MRL should be especially considered since it 
reflects all river stages with only a small degradation in 
Meeeal stabviity and is continuous with MSL, as is discussed 
below. The higher river datums are clearly least stable. 

Regarding continuity, a vertical reference surface, to 
be useful as it passes through the ocean/river transition 
mome, should either be continuous or be constructed from two 
or more intersecting datums so as to avoid any vertical dis- 
Memeinuity. The MRL datum, by virtue of its similar defini- 
miemeana Computation to MSL, 1s continuous through the 
mPansition zone as long as the hourly water levels are 
averaged over the same 19-year period, and indeed becomes 
ome ine MHMRL and MLMRL datums have no counterparts among 
mfeetidal datums; however, MHMRL can be shown to intersect 
MHHW and MHW, and MLMRL intersects MLLW and MLW, thus satisfy- 
ing the continuity requirement. These two river datums 
mmedcieiel to within 0.5 to 0.7 feet at San Francisco but 
never reach MSL in the purely oceanic environment. The MFL 
and MLRL datums also have no counterparts among the tidal 
datums, nor do they intersect any standard tidal planes. 
They appear, therefore, to have limited usefulness in the 


transition zone. MLRL may have some value, however, as a 
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Navigation datum in the Sacramento River during periods of 
very low river level due to its minimum elevation, good 
Stability, and ease of determination. Continuity relation- 


ships between tidal and river datums in the transition zone 
are summarized in Table 6. 

With regard to tidal datum planes, MSL is the only plane 
that is continuous through the transition zone, as previously 
stated. It may be determined in the same manner in both the 
ocean and the river. All other standard tidal planes con- 
verge upon this plane as they are extended upriver. 

From a practical standpoint, the tidal planes become 
more difficult to determine upriver as river stage variations 
become larger. As was indicated earlier, the range of the 
tides decreases with rising river level. At river levels 
greater than about 15 feet at Sacramento, the tides are no 
longer detected reliably as their range diminishes to the 
magnitude of the random noise in the data (0.1 foot). Thus, 
eurins periods of high runoff, individual tides from which 
tide-related datums might otherwise be determined are unde- 
fined. Even before that stage is reached, however, distinc- 
tion between the higher and lower of the high waters, and 
Bisco ot the Low waters, becomes uncertain so that the asso- 
ciated datums, i.e., MHHW and MLLW, cannot be determined. 

MHW and MLW might be determinable under the latter condition, 
however, as suggested in the following paragraph. 

A possible solution to the disappearing tide problem 


that occurs at times of high river stage, and which also 
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eaeteorat, ole Upriver limit of tidal influence, is the crea- 
tion of datum planes (other than MRL) by a separation-addi- 
£i0n procedure. The raw river data can be decomposed into a 
tidal and a non-tidal component using the simple meaning tech- 
nique described earlier in this study. A mean river level 

can be calculated from the raw water-level record (or from 

the non-tidal component as discussed above). From the tides 
that are extracted, a mean tide range can be determined. By 
adding (or subtracting) one-half of the mean tide range to 

the mean river level, for example, an artificial mean high 
water (or mean low water) can be determined. In computing the 
mean tide range, the range of all tides including those that 
have disappeared into the noise background would have to be 
considered in order not to bias the resulting artificial MHW 
level (or MLW level) toward an extreme value. 

A non-technical drawback to this procedure might be un- 
willing acceptance because of the more complex and indirect 
nature of the procedure. Although MHW and MLW could be com- 
puted by this method, MHHW and MLLW cannot because of the dif- 
ferentiation problem noted above unless values of HHW and LLW 
Meam@m@eceansition station are designated according to the 
tide occurrence at the ocean reference station (and tides 


that have disappeared are similarly accounted for). 
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V. SUMMARY 


This study has investigated the problem of extending tid- 
mieeand river datum surfaces through the ocean-river transi- 
tion zone of a navigable river. It presents the results of 
analysis of hourly water-level measurements made over a 19- 
year interval (1959-1977) at two gaging StatLrons lecated on 
the Sacramento River system, one station being the Presidio 
tide gage at San Francisco which represents a purely tidal 
regime and the other at Sacramento representing at times a 
tide/river regime and at other times a purely river regime. 

In order to know how tidal planes may be extended up the 
tidal reach of the Sacramento River, a comparison was made 
of the tide ranges and tide times at Sacramento with San 
Francisco, as a function of river stage at Sacramento. To 
accomplish this, the water-level data at the two stations 
were decomposed essentially into tidal and non-tidal compo- 
nents By use of a 25-hour running mean of the hourly values. 
Wuesmemecedire effectively separated out the prominent semi- 
eeeeneal and diurnal tidal components and left a residual of 
non-tidal river stages, meteorological effects, and presum- 
ably also tidal components of long period but very low ampli- 
tude. Computer programs were developed to determine the 
em@omeceand times of high/low waters to 0.1 foot/0.1 hour 
Mecieesien, and to filter out erroneous data. It was deter-— 


mined that as the river level increases due to increased 
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Bimert, ene ratio of the tide range at Sacramento to that at 


San Francisco decreases, and the travel time of the tide 
wave up the river syeteu increases. At very high river 
Stages, the tides at Sacramento disappear into the back- 
ground river stage noise, which has a magnitude of approxi- 


mately 0.1 foot. Under these conditions the regime at 


Sacramento is purely riverine. 

The six standard tidal datum planes at San Francisco 
Soe eevee MIL, MSL, MEW, AND MLLW) and five river datum 
planes defined and considered in this study (MFL, MHMRL, 
MRL, MLMRL, AND MLRL) were examined for their continuity and 
interrelationships through the ocean-river transition zone. 
iPieaddteton, the stability, or vertical variability, of the 
river planes was determined. Regarding the latter, the low 
water datum planes were found to have the greatest stability. 

Seeere tidal datums considered, only one was found to 
extend continuously through the transition zone -- the MSL 
Sacaum whieh bécomes MRL upriver. This "plane" is easily 
visualized and calculated from digital water-level measure- 
ments made at any location in the transition zone and at the 
Same time has relatively good stability in the riverine en- 
vironment. Usage of this datum would require that MSL and 
MRL be calculated over the same 19-year interval to avoid a 
discontinuity in the plane. 

ferme other tidal datums, MHW and MLW could be carried 


completely through the transition zone by a separation-addi- 


tion procedure involving separating the tidal from the 
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riverine component, determining the mean range of the tidal 
Sempenene (taking into account all of the twice daily tides 
including those ana have merged into the 0.1 foot back- 
ground noise), and adding/subtracting half of this mean 
range (which might be called the mean amplitude) to/from 
the mean of the residual hourly values to obtain a mean high 
water plane or low water plane. These datums effectively 
merge with the MRL datum above aha “seas of the tides. 

The MHHW and MLLW tidal datums could be extended into 
the pure river environment by a similar application of the 
S-pavatton—addition procedure, but this additionally requires 
Maemelfication of which of the two high waters occurring 
daily should be accepted as the higher high water in comput- 
ing the HHW amplitude (a similar argument applies to the LLW 
amplitude). Such an identification procedure 1s necessary 
because when the tidal component is small during higher river 
stages successive values of HHW Cor of LLW) are sometimes 
reversed from the ocean tides recorded at the reference sta- 
tion and at other times cannot be distinguished in recording 
Pema wanplitudes tO a precision of 0.1 foot. The logical 
identification procedure would be to choose those tices that 
are equivalent to the HHW or LLW values determined at the 
Menremence tide Station (San Francisco, in the case considered 
here). Like the MHL/MLW tidal datums, the MHHW/MLLW tidal 
datums also effectively merge upriver with the MRL. 

Concerning the river datums listed above, it has already 


Boepeesseeaoechat MRL is continuous with MSL. The two datums 
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MHMRL and MLMRL have been shown to have no counterpart among 
the tidal datums, although they intersect the MHW and MLW 
planes, eeseetiyely . in the transition zone, and could there- 
mome provide continuity with those planes through the transi- 
tion zone. MFL and MLRL neither are represented in the 

purely tidal regime by a standard tidal plane nor intersect 

a standard tidal plane, hence neither can provide continuity 
fareoueh the river/ocean sae i ellen zone. 

fomeeneeerarad tO the Pacific Coast charting datum, MLLW, 
and the waterfront property boundary datum, MHW, the follow- 
ing summary statements can be made. MLLW can, by a procedure 
of separation-addition and tide identification, be extended 
continuously through the ocean/river transition zone where 
it merges upriver with the MRL datum. MHW can be similarly 
extended continuously through the transition zone by a sepa- 
ration-addition procedure where it also merges upriver wit! 
myeeRo Gatum, or it can be carried through the transition 
zone to intersection with the MHMRL datum thence upriver as 
that datum. There appears to be no way in which MHW can be 
extended as a continuous surface through the transition zone 
SO as to become a low river datum upriver. 

In this study the water-level history at a single gaging 
Seaeronmy 1.c€., Sacramento, was analyzed. It should be recog- 
Mized that in order to carry datum surfaces through the 
ocean/river transition zone with a specified precision at 
any given location, water-level gaging stations would need to 
Pewestabtished at suitable intervals through the transition 


Zio) 
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This study represents an exploratory examination of ver- 
ieeciedatums in the tidal reach of a navigable river. It 
is hoped that ideas and procedures presented here may prove 
useful in future tide and river level analysis and datum 


plane determination. 
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Table 3: RANGE RATIOS BETWEEN SAN FRANCISCO AND 
SACRAMENTO 
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Table 3 (Continued) 
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Table 5: 19-YEAR DATUM PLANE ELEVATIONS AT 
SAN FRANCISCO AND SACRAMENTO 
All datum heights are referenced 


to NGVD 
Sample Mean 95% Confidence Standard 
Datum Plane Heweh t Half-width Deviation 
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tebe S: “DATUM PEANE RELATIONSHIPS 


Ocean Datums 


River Datums MHHW MHW UAE IG: MSL MLW MLLW 


Pe 0 0 0 0 0 0 
MHMRL x x 0 0 0 0 
MRL 0 0 0 Gas 0 0 
MLMRL 0 0 0 0 x K 
Mig LL 0 0 0 0 0 0 
0 = Planes do not intersect in the transition zone 
hMoeenPtanes intersect in the transition zone 

C = Planes are continuous through the transition zone 
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Figure 8: COMPUTER CALCULATED HIGH/LOW WATERS SUPERIMPOSED 


ON RAW DATA CURVES (1-5 OCT 1968) 
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Appendix A: COMPUTER PROGRAMS 


Appendix A-1l 


This computer coding is used to separate the raw 
water-level data into a tidal and a non-tidal component 
By computation of a 25-hour running mean of the hourly 
water levels. "Heights" are the yndividual raw data 
Bemmes. 'Amean'’ 1s the 25-hour running mean value, and 


Sirde” 35 the extracted tidal component. 
RAW WATER-LEVEL DECOMPOSITION PROGRAM 


HH=Height(1)+Height(2)+...+Height(25) 
Do 30 Jsl, K 

Amean(J)=HH/25. 

Tide(J) =Height(J+12)- Amean(J) 

HH=HH - Height (J)+Height(J+25) 


30 Continue 
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Appendix A-2 


ia S computer coding is used to determine when a high 
or low point in the data has been passed. It uses a finite 
difference scheme to give a rough first derivative of the 
curve containing the hourly data points. In this program 
"Tide" 1s the extracted tidal component of the raw water- 
level data. "IC" and "ID" are indicators which determine 


the change in slope of the curve. 


HIGH/LOW WATER INITIAL GUESS PROGRAM 


Do 60 K=il1, LK 
A = Tide(kK) + Tide(K+l1) - Tide(Kt3) - Tide(Kt4) 


PetCo, 2.3 


eG. = Se 
Go to 4 
eee 2c = +) 


ft) scone emue 


If (ID + Ic) 60, 5, 60 


eC = = IC 
ope LDS LC 
- 
+ 
- 


60 Continue 
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Appendix A-3 


The five hourly data points surrounding an initial 
guess at a high or low water level were fitted to a quadra- 


SLC Equation of the form: 


where y is water level and x is time. The normal equations 


for a second-order least-squares fit are: 


Wy 
oO 
Ea 
“be 
ay 
a 
~ 
bh 
+ 
ay 
NO 
M>- 
~ 
fre 
NO 
ti 
M4-= 
eG 
fe 


1=L1 cae - =k 
N N N N 
+ ; a ; Jee aan 
a0 et a el 9 i. ivi 
1=1 1=1 ell izl 
N N N N 
oe e. : | 3 a - yo a 2 
19 at ail xa ? ale 7 ae 14 
ae 1=1 i=l i=1 


For this investigation N always equals 5, and by adjusting 
B@emamdexing the times can always be 1, 2, 3, 4, and 9. 


In this case the normal equations reduce to: 
5 
24.4 P 15a, a 99a, = ) Ne 


15a, a 29485 a 229a,= | XY 


s 
> ee 
294, +225a, E eae, cere) 





Solving these equations simultaneously gives the following 
equations for the coefficients of the quadratic in which 


the water-levels are the only unknowns: 


, 5 5 
= = 2 = 
a. ( S/1L) LCD) oe allele =. y,)/6] 
i=1 ee) 
5 a: 
Z 
= ORs >> Ky. = 55 > XY. )/55)} 
1=1 1=1 
= 5 
= 2 | 
a, = I ) x.y, - ll > y,)/601-(374a,/60) 
1=1 i1=1 
5 
ag = o> 2 3a, - lla, 
1=1 


"Tide" refers to hourly water-levels associated with 
the extracted tidal COmpenents 1.c i ener utdal component 
extracted by use of the 25 hour running mean (given in 
Appendix A-1) applied to either ocean or river water-level 
data. Naa", Nay" and Nag" are the coefficients of the 
quadratic being fitted. "Time" is the differential of the 


quadratic set equal to zero, therefore the time of high/ 


low water. "Height" is the high/low water height. 
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CURVE-FITTING PROGRAM 


Os J + 1 
P24 = Tide(K) + 4.*TideCKtl) + 9.*Tide(K+2) + 16.*Tide(K+3) 


aoremagde (K+ ) 


Meme tade(K) + Tide(K+l) + Tide(K+2) + Tide(K+3) + Tide(Kt4) 
Bie) tade(K) +2.*Tide(Ktl) +3.*Tide(K+2) +4.*Tide(K+t3) 
+5 .*Tide(K+t4) 
wero le) ( (CP 2-11. *FC)/6.)=-( 15.*F2-55."°FL)/ 35. ) 
Bee: 2-11.*FC)/60.-((374./60.)*a2) 


POm= FC/5.=-3.*A1-11*A2 


Time(J) = -Al/(2.*A2) 


Height(J) OMe eshte (0) 9 eleeeh Teme Gul) 22 


Bemec)) = Time(J) + Float(K) 
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Appendix B: GRAPHS OF ANNUAL WATER-LEVELS AT 
SACRAMENTO AND SAN FRANCISCO 


This Appendix consists of graphical printouts of hourly 
water-levels for the 19 water years (1959-1977) studied at 


sacramento, followed by four years of San Francisco data 


(1959-1962) for comparison. Both tidal and non-tidal 


components are shown. 
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Appendix C: MONTHLY AND YEARLY MEAN 
WATER-LEVELS AT SACRAMENTO 

This Appendix contains the monthly and yearly mean 
water-levels for Sacramento derived from hourly heights. 
Given are both the mean of the raw data and the mean of 
the non-tidal component. The closeness of the values 
provides the rationale for eSnipiene river datum planes 
uSing either raw or residual data. The value in paren- 


thesis indicates an incomplete month of data. 
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